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nilEODUCTION 

Let 

(1)  A:k  «  b 

denote  a  system  of  m  linear  equations  In  m  unknowns,  and  let 


(2) 


X  .,  »  Mx  +  k 
ttfl  n 


be  an  Iterative  method  for  finding  the  solution  to  (1).  There  are  two  central 
questions  In  using  such  an  Iterative  method  to  approximate  the  solution  to  a 
set  of  linear  simultaneous  algebraic  equations: 

(I)  Does  the  Iterative  method  (2)  ultimately  converge  to  the  solution  of  (1), 

l.e.  does 

11m  X  =  X  ? 
n 

tt+« 

(II)  What  Is  the  error  e  or  and  error  estimate  E  for  the  difference  between 

n  n 

the  computed  solution  and  the  true  solution,  l.e.  what  Is  such  that 


1 1  ®n  1 1  *  1 1  ^  ^  ^ 


n 


The  first  of  these  two  questions  has  received  considerable  exposure  In 
the  li '.stature.  The  texts  by  Varga  and  Young  give  convergence  criteria 
for  a  variety  of  methods;  moreover,  each  of  these  references  has  an  extensive 
bibliography  on  the  subject. 

Error  estimates,  on  the  othsr  hand,  have  received  much  less  attention.  The 
results  that  have  been  derived  can  be  broken  doim  into  three  categories: 

TYPE  I.  The  error  estimate  E^  assumes  a  knowledge  of  the  inverse  or  an 
approximate  Inverse  of  A. 

TYxE  II.  The  error  estimate  E^  Is  given  only  in  terms  of  previous  iterates, 
for  example  suppose  there  exists  a  computable  constant  a  such 
that  for  some  vector  norm 


n 


x-x  <E  =  a  x  -X  , 
n  ' '  —  n  ' '  n  n-1 


1 


NOLTR  72-189 


TYPE  III.  The  error  Is  given  in  tens  of  two  vectors  t  and  v  which 

n  n 

bound  both  the  solution  x  and  the  n-th  Iterate  x_.  l.e.  at 

n’ 

each  Iteration  u  and  v  are  calculated  such  that  for  aome 
n  n 

partial  ordering,  jc  ,  between  vectors 

u  <  X  <  V  and  u  <  x  <  v^. 
n—  n—  n  n—  —  n 

In  this  paper  we  will  discuss  each  of  these  three  types  of  error  bounds. 

Special  emphasis  will  be  placed  on  the  estimates  of  Type  II  since  the  estimates 

fg) 

of  Type  I  have  recently  been  thoroughly  reviewed  by  Fitzgerald  '  ,  and  the 

estimates  of  Type  III  generally  apply  only  to  very  specialized  pre blame. 

We  will  review  the  articles  that  have  appeared  in  the  literature  which  give 
error  estimates  of  Type  II,  beginning  with  the  historically  Interesting  paper 
of  von  Mlses  and  Pollaczek-Gelrlnger  (17) ,  chronologically  proceeding  to  the 
latest  papers  available,  and  concluding  with  a  numerical  example  using  the 
successive  over-relaxation  (SOR)  method  and  an  error  estimate  of  Albrecht 
NOTATION  AND  BASIC  CONCEPTS 

Throughout  this  paper  we  will  consistently  use  the  terminology  of 
Householder  and  Varga  .  For  completeness  the  definitions  and  concepts 
used  In  this  paper  are  presented  below. 

Let  R  denote  the  real  numbers,  C  the  complex  numbers,  E*"  the  m-dlmenslonal 
vector  space  over  the  field  complex  numbers,  and  GCE*")  the  set  of  all  m  x  m 
complex  matrices.  Then  1 1  •  1 1  :  E®  R  is  a  vector  norm  on  E®  if 

(a)  ! I  X  i 1  ^  0  for  all  X ,  and  | |  x  | |  »  0  if  and  only  If  x  •  0 

(b)  I |X  X  I  I  -  !x|  •  I  I  X  1 1  ,  X  In  C 

(c)  M  X  +  y  I  I  <  I  I  X  1 1  +  I  I  y  1 1 ,  X,  y  in  e“. 

Furthermore  | |  •  | | :  G(E®)  ->  R  is  a  matrix  norm  If,  in  addition  to  (a),  (b),  (c), 

(d)  1 1  AB  1 1  <  1 1  A  1 [  •  1 1  B  1 1  for  all  A  and  B  in  G(E®)  Is 
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satlsf  le  1  ■■  Let 


be  eay  vector  nom  on  i  then  the  equation 


II  All 

x^o  II  X  II 

defines  a  matrix  norm  on  G(E^)  which  is  said  Lo  be  Induced  by  the  vector 
norm  | |  •  | |  on  e“.  A  matrix  norm  { |  •  | |  is  consistent  with  a  vector  norm 

II  •  II  if 

II  II  ill  A  II  •  II  .  II 

for  all  A  in  G(E®,)  and  all  x  in  E™;  moreover  the  matrix  norm  Induced  by  a  vector 

norm  is  consistent  with  that  vector  nor.'.  Let  p(A)  denote  the  spectral  radius  of 

A,  then  pCa)  ^  1 1  A  1 1  for  any  matrix  norm. 

*  * 

Now  let  X  £ad  A  respectively  denote  the  conjugate  transpose  of  the  vector 
X  and  the  .atrix  A,  where  A  =  and  x  »  •  Three  widely  used  vector  nonns 

and  their  induced  matrix  norms  are: 


m 

m 

I  1  "<4  1 

,  II  A  1 L  =  max  ^ 

1=1  ^ 

l<J<m  i=l 

J* 

\X  X  = 

^ 1  1  1^.  1 !  A 

91=1 

m 

max  1  Xj^ 

1  .  1 1  A  1 1^  =  max  1 

/p(A*  A) 


l<l<m 


j=l 


It  is  well  known  (Varga  that  the  sequence  x^  from  (2)  converges 

for  any  x^  if  and  only  if  p(M)  <  I  If  A  is  nonsingule:  ..nd  the  method  (2)  is 
derived  from  the  system  (1),  i.e,  x  =  A  Ic  the  only  solution  to  x  »  Mx  +  k, 
then  X  =  lim  x 

n 

n 

is  the  unique  solution  to  (1). 
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The  three  most  commonly  used  Iterative  methods  are  the  Jacobi*,  Gausa-Seidel**, 
and  successive  over-relaxation  methods.  Each  of  these  three  methods  will  now 
be  expressed  in  the  form  (2).  Decompose  A  into  a  diagonal  matrix  D,  a  lower 
triangular  matrix  E,  and  an  upper  triangular  matrix  F  such  that 
(3)  A  -  D  -  E  -  F. 
then  the  Jacobi  method 

-  »j,(»)  *  ‘i 

can  be  written 

vi  ■  ^  ^Vi  + 

so  that  the  Jacobi  matrix  B  is  D  ^(E  +  F).  The  Gauss-Seldel  method 


(5) 


*l,(tri-l) 


1-1 

jll 


ji+i  "i  (") 


can  be  written 

x^j^  -  (D  -  E)“^Fx^  +  (D  -  E)"^b 

so  that  the  Gauss-Seldel  matrix  is  (D  -  E)  ^F.  The  successive  over-relaxation 


method 


(6)  X, 


i-1 

-  I 


-  I 


i,(n+l)  a^^  ij  J,(n+1)  ij  J  ,  (n)  i 


+  b.  -  ('i)-l)x. 


can  be  written 

X  ^  -  (D  -  a'E)'^  [(1  -  w)D  +  (oF]  x  +  a)(D  -  u)E)~^b 
nri  n 

ejacobl  method  is  also  called  the  point  total  step  method,  the 
method  of  simultaneous  displacements,  and  the  Richardson  Iterative  method. 
Note  the  total  step  method  is  a  translation  of  the  German  word 
Gesamtschrittverfahren. 


**The  Gauss-Seidel  method  is  also  called  the  point  single  step  method, 
the  method  of  successive  displacements,  and  the  Llebmann  method.  Note  that 
single  step  method  is  a  translation  of  the  German  word  Einzelshrlttverfahren. 
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so  that  the  successive  over-relaxation  matrix  is  (D  -  uE)"^  [(1  -  (d)D  +  wP] . 

We  note  that  for  (o  *  1  the  successive  over-relaxation  method  reduces  to  the  Gauss- 
Seidel  method.  Young  shows  that  under  certain  conditions  there  exists  an 

optimum  over-relaxation  parameter,  1  <  w  <  2,  such  that  p(L  )  <1,  aud  that  the 

(i) 

successive  over-relaxation  method  then  converges  faster  than  either  the  Jacobi  or 
Gauss-Seidel  methods. 

We  now  proceed  to  discuss  the  error  estimates  of  types  I,  II,  ar:*  III. 

ERROR  ESTIMATES  OF  TYPE  I 

(8) 

Recently  Fitzgerald  '  •'  developed  bounds  for  the  error  in  a  computed 
Inverse  of  a  matrix  and  for  the  approximate  solution  of  Ax  ■  b.  These  bounds  are 
of  the  first  type  listed  above  in  that  he  assumes  the  availability  of  an  approxi¬ 
mate  Inverse  X  of  A.  Methods  in  which  these  bounds  would  be  most  applicable 
compute  an  approximate  inverse  X  of  A  and  let  y  ■  Xb.  The  residual  I  -  XA  or 
I  -  AX  is  then  computed  and  from  these  a  bound  on  the  error  e  ■  x  -  y  is  found 
using  the  well  known  Inequality 

i  1  V  -  X  I  I  <  I  I  I  I  •  I  I  Ay  -  b  I  I 

where 

II  A-i  II  < — UjlU - 

l-ll  I-AXII 

if  I  I  I  -  AX  I  I  <  1. 

If  this  error  is  too  large  X  may  be  improved  until  the  error  is  within  acceptable 
bounds.  When  one  uses  such  a  method  to  solve  for  the  solution  vec^-or  x  of  AX  ■  b 
for  just  one  b,  he  is  actually  computing  much  more  than  he  neSds,  namely  A 

Fitzgerald  feels  that  in  general  it  is  futile  to  expect  to  find  a  good  error 
bound  without  some  knowledge  of  the  inve“se  of  A;  however,  one  would  hope  that  if 
A  had  a  "special  enough"  structure  it  would  be  possible  to  find  an  error  estimate 
without  knowing  anything  about  A  that  is  find  error  bounds  of  the  second 
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or  third  types  above.  As  we  shall  see,  this  is  true  for  some  frequently  used 
Methods  eaplcqred  on  classes  of  matrices  which  often  arise  in  the  solution  of 
problems  derived  from  physical  phenomena. 

ERROR  ESTDMIES  OF  TYPE  II 

He  now  present  the  error  estimates  of  the  form 


X  -  X 


<  O  X  -  X  - 

' '  n  n-1 


Although  some  of  these  bounds  were  originally  neither  presented  nor  proved  in  this 
form,  for  siiiq>llclty  and  self-consistency  they  will  be  given  using  the  results  of 

/1Q\ 

the  following  theorem  due  to  Welsslnger  . 

Theorem  I.  Let  1 1  •  j  |  be  a  matrix  norm  consistent  with  the  vector  norm 
II  •  II,  and  let  x^^ =  Mx^  +  k  be  a  method  derived  from  the  nonsingular 
system  Ax  »  b.  If  | |  M  | |  <  1  then 


X  =  11m  X  ,  and 


o>  II  X  -  II  < 


H 


X  -  X  , 

n  n-1 


1  -  II  M  II 

Proof:  Since  p(M)  <  ||  H  ||  <1  the  method  converges  to  the  unique  solution 

of  AX  *  b.  For  the  proof  of  the  second  statement  let 

e  ■  X  -  X  ,  and 
n  n* 

n  n  n-1 

Then  since  6  , .  ■  M5  and  p(M)  <  1,  it  can  be  easily  shown  that 
n+i  n 

e_  -  (I  - 


and  since 


is  a  consistent  matrix  norm 

.-1 


I  X  -  xjl  <  II  (I  -  Mr  II  •  II  M  II  .  II  I! 
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••1  *  dL 

Now  again  since  p(M)  <1  we  can  expand  (I  -  M)  as  I  M  ,  then  using 


1=0 


(I  -  M)"^  M  <  f  I !  M  M  ^  = _ i_ 

1=0 


1  -  M  M  II 

we  have  the  result  of  the  theorem. 

We  can  see  from  Theorem  I  that  If  we  can  find  an  easily  calculable,  consis¬ 
tent  matrix  norm  such  that  | |  M  | |  <1,  then  we  have  a  computable  error  bound. 
This  Is  In  fact  what  has  been  done  for  the  Jacobi,  Gauss-Seldel,  and  successive 
ever-relaxatlon  methods  using  the  three  norms  ||  •  ||j^,  ||  *  1 1 2»  II  *  Moo* 
We  must  note  however  that  In  general  any  or  all  of  the  above  norms  may  give 
II  M  II  >1,  precluding  the  use  of  the  error  bound  (7),  but  the  method  may  still 
converge. 

We  now  present  the  error  bounds  of  Type  II  that  have  been  derived.  In  all 

that  follows  we  assume  we  are  solving  the  m-th  order  system  Ax  =  b,  where 

A  «  ^  “  (b^) »  and  x  =  (x^) ,  by  the  method  “  ^n  ^  where 

X  ■  (x^  .  >)  and  M  is  matrix  of  order  m  which  may  be  one  of  the  matrices  B,  L.  , 
n  i,(n)  ^ 

or  L  corresponding  to  the  Jacobi,  Gauss-Seldel,  or  successive  over-relaxation 


methods  respectively. 


.(17) 


In  1929  von  Mises  and  Pollaczek-Geiringer '  proved  that  if 

<  V  <  1  for  j  =  1 ,  ...  ,  m 


m 

r 

1*1 


!ll 


*11 


then  an  error  estimate  for  the  Jacobi  method  is 


m 

I 

1=1 


'‘i,(n+l)  "  ’‘i 


ra 

<  y  I 
1-y  1=1 


^i,(n+l)  ■  i,(n) 
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Equivalently  using  the  terminology  and  the  result  of  Theorem  I  we  have: 
Ify  =  j|Bl|j^<l,  then 


^n+l  ■  ^ 

*n+l  ~  * 

1  1-P 

1' 

No  other  results  were  obtained  until  1942  when  Collatz  showed  that  for 
the  Jacobi  method,  if 

m, 

af 

— ^  5  <  1,  for  i=l,  ...  ,  m 


I 

j=l 


ii 


then  Max 

1  <  i  <  m  I  X 


i,(n+l)  ^i 


'lax 

<  1  <  1  <  m  j  X 

l-€ 


i,{n+l)  ■  *l,(n) 


He  also  noted  the  same  bound  held  for  the  Gauss-Seldel  method.  Again  we  restate 
this  bound  in  the  form  of  Theorem  I  to  get:  If  5  =  I j  B  | [  <1  then 


1  ^n+1  "  ^  1 

<  j  j  ^n+1  "  ’^n 

00  1-r  1  I 

Again  there  was  a  long  lapse  until  1951  when  Sassenfeld presented  two 
results  for  the  Gauss-Seidel  method,  naraelv. 

Criterion  I.  Let 


1  i-l 

a.  =  -  T 

^  ^ii  j=l 


‘ij 


m 


cx .  4-  I 

J  j=l+] 


,  i=l. 


m 


and  a  =  Max  a,.  If  a<l  then 


Max  I  X . 


i,(n+l)  i 


-  X .  I  <  a  Max 


1-a  i 


^i,(a+l)  ^i,(n) 


Criterion  II.  Let 

=  — — 


i-1 

r 

Lj=i 


'Max  Q  \  ,  P 


,  i=l  ,  ...  ,  m 
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and  S  =  MaxB..  If  B  <  1,  then 
i  ^ 


Max  X .  ,  .  -  V  -  X .  <  r— r-  Max  x .  ,  . ,  v  -  x ,  ,  ^ 

I  i,(n+l)  1  !  1-B  ^  l,(n) 


We  observe  that  Criterion  I  is  the  stronger  of  the  two  since  a<B.  Both  of 

Sassenfeld's  criteria  will  be  proven  as  corollaries  to  the  following  more  general 

theorem  on  the  successive  over-rel;ixation  method  which  uses  Theorem  I  with  an 

explicit  representation  of  and  the  norm  j |  •  | |^.  Referring  to  (6)  we  can 

write  the  successive  over-relaxacion  matrix  I  =  (s..)  as 

tu  13 

s^j^  =  -(u  -  1) 

s,.  =  u)  a,. 

IJ  Ijj  j  2,  ...  ,  m 

and  for  1=2,  ...  ,  m 

i  -i 
^  k=l 


®ii  = 


if  T 

=  -  - —  (u  -  l)a.  .  +  u)  I  a.,  s,  .  I 
^ii  [  k=l 


"ij  =  -  (a.j  +  u,  E  a^^s^.),  j=i+l. 


Theorem  II.  Let  Y  =  I  and  y  =  Max  Y.,then  if  y<1.  an  error 


estimate  for  the  successive  over-relaxation  method  is 


i  i  ''i,(n+l)  '  ""i  '  i  '  ’'i,(n+l)  "''i,(n) 


Prcof:  We  note  that  y  =  I  L  |  .  Then  If  <l  we  have 

U  CO 


n+1 


X  , ,  -  X  I 

n+1  n 


which  is  the  conclus-oon  of  the  theorem. 
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Let 

I  u  I 

*11  J-1 


r  1-1 

_1 _  U  I  I  I  +  I  ”  l)aj 

*u  I 

“I  I  a^^  1  .  1*2,  ...  ,  n 

j-1+1  J 


and  5  *  Max  5,. 

1  1 

If  5  <  1  then  for  the  successive  over-relaxation  method,  ve  have 

I  ’‘l.(nH)  -  *1  I  ^  I- 

Proof:  We  need  only  to  show  that  y  <  6.  For  1  >  1  we  can  write 


ri  -  !  I  Sjj 

j-1 


‘  1-1  1-1 

u  /  y  a.i  s,  . 

jii  w  “‘‘J 


+  (.  -  1).  +  «  I  a 

k«]  '* 


+  u  I 


m  i-1 

I  a  +  0)  y  a  s 
J*l+1  k-1 


I  L 

^  [“  ' 

I  L 


I  I  I  s. ,  I  +  I  (lo  -  Da..  I  +  0)  y  |a  J 


Ik  '  ^  '  kj 

j-1  *' 


j-1+1 


I  ay  I  Tj  +  I  («  -  Day  I  +  a  ^  I  ay  I 


Now  ■  Yj^*  hence  1  5^^  and  y  <  6. 
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CorollTy  II.  Let 

V,  •  6, 


‘  *  I  »L  I  f* 


I  +  I  (»  -  1)  I 


+  «  I  I  »ii  I  • 

J-1+1 


and  V  »  Max  v\ 


If  V  <  1  then  for  the  successive  over'-relaxation  aethod,  «e  have 

Ms*  I  \  ~  I  <  M«*  I  \  \  I  • 

^  '  i,(iH-l)  i  '  -  1-v  i  '  l.(iH-l)  i,(n)  ' 

Proof:  This  corollary  obviously  folloirs  froa  Corollary  I  since  v  <  6. 

Corollary  III.  Criterion  1  and  II  of  Sassenfeld. 

Proof:  Let  u  «  1  In  Corollaries  I  and  II. 

Dueck  presented  the  following  error  eatlaate  for  the  Gauss^Seidel  aetbod 
which  is  slightly  better  than  Collatz's  estlaate  but  not  as  good  as  Sassenfeld 's> 
Let  A2  be  the  upper  triangular  part  of  the  Jacobi  matrix  B.  If  1 1  B  1 <  1 
then 


II  A,  I! 


“  * 

n+1 


l-ll  B|i 


ttfl  n 


Let  be  the  lower  triangular  part  of  B,  then  Dueck  proved  this  estlaate  by 

noting  that  an  equivalent  fomnilation  of  the  Gauss-Seldel  method  Is 

X  . ,  -  A-x  ,  ■  A-x^  +  k 
n+1  1  r+l  2  n 
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X  -  X 


"  Vl^  ^  "  V 


(I  -  A^)(x  -  x^j)  »  k^ix  -  x^) 


“  ^2<*  -  W  ^2^Vl  ‘  V 


(I  -  Aj  -  AjXx  -  x^^)  =  k^(x^j_  -  x^) 


*  -  Vi  “  ^2^Vi  -  V* 


Therefore  by  Theorem  I 


*2  II. 


*n+1  ^  I-  - - I  1  *..^1  “  * 

iH-1  -  1  _  11  B  II  •  n+1  n 


Feldnan  ‘'found  an  error  estimate  for  the  Gauss-Seldel  method  mhich  is  comparable 
to  that  of  Dueck.  As  before  let  A^  be  the  lower  triangular  part  and  let  A2  be 


the  upper  triangular  part  of  B.  Let 


1-  !i  Aj  II. 


l-ilAjL  " 

I  mil  A2  11,  if  II  Aj  |!„  =  1. 


il  A2  !Lif  II  A^  IL#1 


If  1(1  <  1  an  error  estimate  for  the  Gauss-Seidel  method  is 


n+1  n  '  ' «  • '  n+1  n 


Following  the  proof  given  by  Feldman,  we  have 


*  (I  -  A^)  "A2 


M  h  IL  ^  II  (I  -  A^)'^  !L  I 


'2  "«• 
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is  the  Jacobi  method  for  the  solution  of 

(8)  (I  -  T)y  »  a 

_i, 

and  the  solution  of  Ax  -  b  is  x  =  D  ^y.  However  T  like  A  is  Hemltlan  and  posi¬ 
tive  definite  while  the  Jacobi  matrix  B  is  generally  not.  Using  the  above 
notation,  Albrecht's  error  estimate  for  the  solution  to  (8)  by  the  successive 
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over-relaxation  aethod  when  A  la  Heraitlan*  podtlve  definite,  and  2-Qrcllc  la 

II  W''  lli^*  II  Vi-’'.  II2 

idiere  ■. 

,  .,,Vl|l|l3(»^*l|7|l^*V0-||T|l3*Vl|l|l^to^*l|T||^ 

and 

y  »  1  -  —  and  p  - - 1. 

u  u 

Soir  since  x  «  D  y,  D  *  diagonal  (d^),  an  error  estlaate  for  the  aolntlon  of 
Ac  *  b  is 

»>  II  Vi-'lla-^^^ll  Vi-'iiz 

The  only  thing  rasainlng  is  to  calculate  I  {  T  \\^,  hut  T  is  Hemitian  hence 
j|  T  II2  *  P(T)*  Varga^^^*  ^^glves  nethods  for  finding  <■  (T)  iihich  are 

based  on  the  iterative  srhene  (2)  and  are  thus  ^uite  easy  to  inpleaent  since  the 
saae  scheae  is  used  *:o  solve  the  set  of  equations. 

Young shows  t'lat  in  the  case  of  2-cyclic  aatrices  (i.e.  aatrices  with 
Young's  Property  A)  if  the  spectral  radius  p  of  the  Jacobi  aatrix  is  less  than 
one  then  the  optinua  over-relaxation  factor  is  related  to  p  by 


2 


and  c(L)  »  u  -  1, 

where  L  is  the  associatrJ  over-relaxation  matrix.  Hence  if  p  1  we  have 

0  <  u  <  2  and  p(L)  <  1.  Now  in  our  case  where  A  is  Hermitian,  positive  definite, 

and  2-cyclic,  the  over-relaxation  matrix  P  for  the  system  (8)  is  similar  to  the 

matrix  L  (see  (6))  for  the  system  (1)  by 
u 

P  =  d\  D~’' 

(i)  u) 
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T  - 

therefore  p(P^}  »  p(L^)  end  p(T)  »  p(B). 

Slallerly  denote  the  consteet  vector  fc  of  (2)  hy  fc^  for  and  by  for 
If  as  before  ne  let 

y.  - 


fc  -  D 
P 


and 

then 


t  ,=L»  +fc,* 
3+J  •  n  X 


y--»Py  +fc. 

■'nrx  arn  p 


y  »  dSs  . 
n  n 


Therefore  by  (9)  an  error  estinate  for  the  nethod  «Lz  -flris 

nt-l  an  i 


(10) 


nt-1 


-  I  il,  s  xyij»  ^  II  -  o  II,. 


■n<-l 


or  a  less  desirable  estiaate  is 

(11) 


II  Vl  -  *  liz  ‘  V  II 


‘n¥l  n  "2 


II. 


To  use  Albrecht's  estlaate  «e  have  two  alternatives,  lie  can  either  solve 
the  systea  (8)  for  y  using  the  error  estiaate  (9),  ard  after  sufficient  conver- 
gence  let  z  »  D  y,  or  «e  can  solve  the  original  systea  (1),  using  the  error 
estlaate  (1).  The  two  nethods  are  equivalent  and  neglecting  any  computational 
aspects  such  as  round-off  error  should  give  the  sane  results.  In  the  *zanple 
at  the  end  of  this  paper  we  decided  to  take  the  second  approach  only  for  the 
reasoi  that  for  our  case  it  was  ecsier  to  program. 

(19) 

We  now  turn  to  a  rather  singular  result  by  Weinberger  which  gives  an 

error  bound  tihose  range  is  determined  from  a  maximum  principle.  Specifically  if 

the  real  matrix  M  of  the  iterative  method  x  *  Mx  +  fc  is  synmetrlc  let 

n+1  n  ^ 
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■xn  7^-lM 


*«•*.- Vi« 


*- V 


-  II  Vi  Hi- 


•  *  Cl*.- 

««4  If  •  II  llj. 

then  If  pOO  <  1  "  e  <  1  the  tmage  of  pocslhle  values  of  f  |  e 


il 


2 

2 


is  equal  to  the  raage  of 

+ _ CS*  _  /  [rt>^-28#  +  Yl 

(1-v)^  ((o- B)#  -  (B  -  y)1^  ' 

on  the  interval 


(1  -  g)B  4-  Y  <*<!_£ 

(1  -  e)«  +  B  i  ^  i  1  -  c* 

The  reqniranent  that  M  be  sYHuetric  generally  linits  this  result  to  the  solution 
of  re^tl  linear  systens  bj  the  Jacobi  nethod 


'nfl 


derived  fron  the  transfomed  systen  (8). 

tl3) 

Finally  Schroeder  used  a  theoren  concemiqg  an  abstract  Iteration  process 

u  *  Zu  to  derive  error  bounds  for  the  Jacobi  and  Gauss-Seidel  nethods  which 
nfl  n 

are  essentially  the  same  as  the  bounds  of  Collatz  acd  Sassenfeld  respectively. 
COmEHTS  ON  THE  ERROR  ESTDfBTES  OF  TYPE  II 

All  of  the  estinates  except  Albrecht's  and  possibly  Weinberger's  only  hold 
if  1 1  M  1 1  <  1  where  1 1  *  1 1  is  sone  consistent  matrix  nora  and  M  is  the  aatrix 
from  the  iterative  process  ■*"  k.  However  the  condition  for  convergence 

is  p(M)  <1,  whereas  it  is  possible  to  have  p(M)  <  1  <  1 1  M  1 1 .  We  would  then 
have  a  convergent  method  but  no  usable  error  estimate.  In  particular  for  the 
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fcowadi  of  voo  !flscs  and  fioll«cxck-€elriqier,  Collacx,  SassenfeU,  the  eutiior** 
Tbcoroi  II,  and  Dneck,  it  is  aecessazy  that  the  aatrlx  A  of  ha  «  h  be  strictly 

dlaagaally  doatoaBt.  i.e. 

a 

\  I  *jj  I  I  I  »  1*1» 
jJl 

This  is  onfortimate  slace  the  aatriccs  that  are  generated  by  aaoy  physical 
problcas,  e.g.  the  nuaericsl  solution  of  elliptic  differential  equations,  do  not 
enjoy  thiii  ptcperty. 

On  the  other  hand  Albrecht's  estiaate  for  the  successive  over-relaxation 
aethcd  relies  only  on  the  spectral  radius  of  the  aatrix  M,  but  it  requires  that 
A  be  positive  definite,  Beraitian,  and  2-cyclic.  In  the  example  at  the  end  of 
this  paper  ue  solve  a  set  of  linear  equations  derived  froa  the  nuaerical  solution 
of  Poisson's  equation  by  the  successive  over-relaxation  aethod.  For  this  exanple 
t*»e  only  estiaate  that  is  applicable  is  Albrecht's. 

ERBOK  ESTIMATES  OF  TTPE  III 

Schrocder^^^^and  Albrecht^^^'^^^'^^^give  error  estiaates  of  Type  III  foi 
nonotone  iterative  nethods.  A  description  of  either  aonotone  aethods  theaselves 
or  conditions  under  which  such  a  aethoa  will  produce  a  aonotone  sequence  of 
vectors  which  con/erge  to  the  solution  of  the  system  of  linear  equations  is 

(14) 

beyond  the  scope  of  this  paper.  However,  Schroeder  gives  a  derivation  of 
aonotone  methods  as  well  as  sufficient  conditions  for  the  convergence  of  these 
methods . 

Given  a  monotone  iterative  method,  it  is  usually  easy  to  write  error  bounds 

for  this  method.  For  example  let  <  be  the  component -wise  partial  ordering 

between  vectors,  x  the  solution  vector  of  the  set  of  equations,  and  x^  and  y^  be 

two  vectors  such  that  x  1  x  <  y  ,  then  if  the  method  gives  successive  iterates 

o  o 
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and  7^  such  that 

*c  -  *1  ^  ^  ^  ^ 

then  an  error  eatlnate  for  z-x  orj-y  la 

n  n 

I  -l  -  ^.(.O  '  ^  •  ^.(«)  '  'i.M  I  •  *“  . . 

Thla  exaaple  la  not  neant  to  give  the  noat  general  eatlaatea  for  all  nonotone 
nethoda.  It  la  hooever  repreaentatlve  of  the  type  of  error  eatlnatea  obtainable 
with  theae  nethoda. 

There  are  a  fen  dlatlnct  dlaadvantagea  to  nonotone  Iterative  methods.  Flrat* 
the  range  of  applicability  of  theae  nethoda  la  amall;  aecond*  it  la  usually  diffi¬ 
cult  to  find  initial  valuea  with  the  desired  properties;  and  third,  the  iteration 
itself  is  more  coaqillcated,  often  as  In  the  case  of  our  example  requiring  tvo  or 
nore  sequences  which  converge  to  the  solution.  On  the  other  hand,  termination 
criteria  for  such  iterations  are  easily  determined  as  the  Iterates  give  both 
upper  and  lower  bounds  on  the  solution. 

NUMERICAL  EXAMPLE— SOLUTION  OF  THE  POISSON  EQUATION 

In  this  section  we  present  an  example  of  Albrecht's  error  estimate  for  the 
solution  of  the  set  of  linear  equations  derived  from  the  discretized  Poisson's 
equation.  The  error  estimate  was  programmed  into  the  computer  code  MACNOL 
(Marker  and  Cell  Method  of  the  Naval  Ordnance  Laboratory)  which  solves  Incom¬ 
pressible,  viscous,  Intilal  value,  fluid  flow  problems  by  the  marker  and  cell 
method.  The  MACNOL  code  is  a  modification  of  the  MACYL  code  of  Prltchett^^^^ 

One  modification  of  MACYL  which  was  incorporated  into  MACNOL  was  to  reprogram  the 
routine  which  solves  the  discretized  Poisson's  equation  to  use  the  successive 
over-relaxation  method  Instead  of  the  Gauss-Seldel  method. 
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As  -j.  note  of  Interest  we  would  like  to  report  that  for  the  sample  problem 
we  solved,  the  successive  over-relaxation  method  converged  about  100  times  faster 
than  the  Gauss-Seldel  method, 

A  brief  derivation  of  the  set  of  linear  equations  arising  from  the  finite 
difference  solution  of  the  Poisson  equation  will  now  be  given.  Using  a  cylindri¬ 
cal  coordinate  system  let  r  be  the  radial  dimension  and  z  be  the  vertical 
dimension.  We  denote  a  finite  difference  cell  by  the  indices  1  and  j  where  1 
varies  with  radius  and  j  varies  with  height. 

Figure  1  Illustrates  this  nomenclature: 

Ar^  "  radial  dimension  of  cell  i,j 


vertical  dimension  of  cell  l,j 

distance  from  axis  to  center  of  cell  i,j 

distance  from  bottom  of  mesh  to  center  of  cell  i,j 
distance  from  axis  to  inner  boundary  of  cell  i,j 

distance  from  axis  to  outer  boundary  of  cell  ItJ 

distance  from  bottom  of  mesh  to  lower  boundary  of  cell  l,j 

distance  from  bottom  of  mesh  to  upper  boundary  of  cell  i,j 


“  ’'j+1  ■  "“j 

Figure  2  illustrates  the  computing  mesh.  Rows  1  and  N  and  columns  1  and  M 
are  physically  fictitious.  They  are  used  for  convenience  in  representing  the 
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FIG.  2  THE  COMPUTING  MESH 
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governing  finite  difference  equations.  The  dooain  of  the  problem  need  not  be 
defined  on  each  cell  of  the  mesh;  by  suitable  flagging  we  can  represent  free 
surfaces.  Figure  3  gives  an  example  of  such  cell  flagging.  The  ccmiputation  is 
only  done  in  full  (F)  cells;  the  boundary  conditions  are  applied  on  the  surface 
(S)  cells  as  well  as  the  axis,  wall,  floor,  and  ceiling.  The  empty  (E)  cells  are 
only  used  ■*’1  the  calculation  to  specify  boundary  conditions.  Let  ^  denote  the 
ratio  of  pressure  to  constant  density  and  let  g  be  the  acceleration  due  to 
gravity.  We  wish  to  solve  the  Poisson  equation  in  cylindrical  coordinates 


^  2  ^  r  3r  ^  .  2 
3r  3z 


with  the  mixed  Neumann-Dl rich let  boundary  conditions 


^  =0  on  the  axis 
3r 


=  w(z)  on  the  wall 
0  r 


ii. 

3z 


=  g  on  the  floor  and  ceiling 


and  is  prescribed  on  the  surface, 

where  R(r,z)  and  w(z)  are  known.  Let  ij)  be  defined  at  cell  centers.  The  MACNOL 
code  solves  the  following  discretized  form  of  the  above  problem  for  (J) 


(12) 


<*l,j  -  ‘*i.j  -  *1-1, j>] 
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FIG.  3  EXAMPLE  OF  CELL  FLAGGING 
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with  the  boondazy  conditions 


*1.1  '  *2.1' 


*M,1  “  *H-1,J  ■* 


*1.1  '  *1,2  -  ***3/2 


♦l,»  •  ♦l,*-!  ♦  »‘Vl/2 


and  4.  .  Is  given  if  i,j  is  a  surface  cell,  and  where  R  ,  and  w.  are  regarded 

i»j  j 

as  known. 

Hhere  a|ypllcable  substitute  the  boundary  conditions  (13)  into  the  equations 
(12).  The  result  is  a  set  of  linear  equations  in  the  unknowns  ,»  where  l,j 
ranges  over  the  set  of  full  cells.  Denote  this  set  of  equations  by  DPE 
(Discretized  Poisson's  Equation).  If  the  set  of  equations  DPE  is  naturally 
ordered,  i.e.  ordered  by  the  rows  or  colunns  of  the  finite  difference  •f»h,  then 
its  aatriz  A  is  symnetric,  positive  definite,  and  2-cycllc.  The  oatrlz  A  is  not, 
however,  strictly  diagonally  dominant,  although  it  is  irreducible  and 

m 

I  I  I  ;  Uul.  i-i.”-.  • 

jJi 

with  inequality  for  at  least  one  1.  That  is,  A  is  irreducibly  diagonally 
dominant . 

To  check  the  usefulness  of  Albrecht's  estimates  for  use  in  the  HACNOL 
code,  ve  programmed  four  similar  sample  problems  whose  solutions  are  known 
exactly.  The  four  problems  are  broken  down  into  two  sets  of  two  problems 
each.  The  first  problem  in  each  set  uses  an  unevenly  spaced  mesh,  the 
second  uses  an  evenly  spaced  mesh.  The  first  set  of  problems  is  simply 
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water  etendlng  etiU  la  a  ^Uodrlcal  toA.  For  the  eecoad  oet  of  ymUear,  la 
order  to  test  the  error  estloeteo  on  s  non-rectsmalsr  oeA*  ve  iotrodacod  the 
artificial  problca  of  an  axl-oynnetrlc  nndeinnter  cavity  nhooe  surface  nos  at 
hydrostatic  pressure.  The  solntion  to  all  of  these  problens  Is  sinply  hydrostatic 


pressure. 

To  reiterate,  Problen  1  uses  an  unevenly  spaced  rectaniular  nsah;  Problen 
2  .-ses  an  evenly  spaced  rectangular  nesh;  Problen  3  uses  an  uneveidy  spaced 
non-rectangular  nesh;  and  Problen  4  uses  an  evenly  spaced  non-rectangular  nesh. 

All  four  problena  use  a  nesh  of  dlneasion  56  In  the  radial  dlrectlm  and 
122  In  the  vertical  direction.  When  the  nesh  is  eveidy  spaced 

Ar  *  Az  •  .5. 

For  the  unevenly  spaced  problens  Ar^  and  Az^  are  given  In  Table  1.  Figure  4 
shows  the  cell  flagging  for  Problens  1  and  2,  and  Figure  5  shows  the  cell 
flagging  for  Problens  3  and  4. 

The  error  estlnates  given  are  (10)  and  (11)  of  Albrecht.  Let  D  be  the 


diagonal  of  DPE  then  fron  (10)  we  have 


II  -♦  11,1  a  II  (*  ,  -  *  )  II, 


where 


A  II  II. 


and  from  (11) 
(15) 

where 


il  Vl-*  lUlB  II 
8.J  II  d-1'2  IIj  II  llj 


We  also  investigated  the  possibility  of  using 


o(L  ) 

(16) 


1  -  P(L  ) 

(i) 


-  *Jl9  -  «  -  4 


^ttfl  ""n  "2 


n+1  -"n  "  2 
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TABLE  I  DEFINITION  OF  MESH  FOB  FROBLEMS  1  AND  3 


Nail73>1t» 


000000 000000000000000  000 


^  S*  000  0000^^000000  0  0  00  0^00  000 

0000000 00 00000000 0000 0000000i 


^•0<'00»00-0 00  000-0000’00  000000 


—  cccccrccrccrcc  00  «r  e^cc4rcciirc«’^«ecc« 


E  «  c  c^c  «grir'<c  «w 
r#fr  0000  00^m^ 


99-0  0»€  ^0  0  00000  00^0  0000^9  #  C  <» 
•  00000  0  0  00  0  00  00000000  00  000 


T  0-€0*<0  <0  cr'Cr  crcrcrcrcccccccc  c>c  i£-c  c  c  c  C'C«c  •€  ^c-c  'Oc  C'C«c  c  c  cc  c  c«C'C'C  c 

•  *  ^^^00  0  00000  0^^  mm  00  0^90*0000-00  0^90*^090  ^00  00  0000*0900  0909 

mm  0  0  0009  00190909  09 


rrr*‘rrirrrirrrtf*rrrtf'#rr««'c^ccccccceeccccc«eccccccc  er  c  c  c  c  e  c 

**********  *  ••*••  ^  BLitfirrrirc  crrccccecc 

^  ^9  90  0  0  000  0 
—  ^999 
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t4 


II 


9  •  •  •  !•  »•  tt  >•  M  *•  V  )«  H  ia  *•  •»  «•  4ft  «•  ft*  ft4 

FIG.  4  CELL  FLAGGING  FOR  PROBLEMS  1  AND  2 
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0 

0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 


0  0 
0  0 
0  0 
0  0 
0  0 
0  0 


0  0  t 
0  0  9 
0  0  9 
0  0  I 


^000900^ 


*  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  •••  9  ^  • 

^099  0  000  0  00  00  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0 

999900000000000000000000000000* 
^99999000000000009***00000*0*^** 
999990000000000000*00000000*00^** 
9999999m000000000000*00*0000***^** 
^999999990900000000000000000000000* 
999999009090900000000000*00000*^^*^** 
999990m90990000*0******»****************** 
9999999090*************************’*******^*^ 
^  ^  m  ^  •  m  m  9  9  9  9  m  m  0  9  0  0  0  0  0  *  *  9  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * 

I  9  9  9  9  9  9  9  9  9  0  0  0  0  9  0  9  0  0  *  •  0  0  *******  0  ****************  * 

^9990999999000000000000*00***0****************** 

9999999999900000000*0************************** 
99999999999  0  0000m00*********0******************* 
^99999999999999099*90*************************** 
999999999990009009000*************************** 
9999999999900090000^0*************************** 
999999999990000000**0*0************************* 
99999999999000900000**********0***************** 
9999999900909990000*00************************** 

9  9  9  9  9  9  9  9  0  0  9  0  0  0  0  0  9  *  0  ***********************  *^*^*^*^*^* 
99999999990999009999*9************************** 
99999m999990099999**0*************************** 
999m»9999999999909999**9***0******************* 
9  9  9  ,  9  9  r  9  9  9  9  0  9  9*99******************************* 

999999W0999«****************************tlll9ll* 

999999999990009000999************************** 
9  9  9  9  9  9  9  9  9  9  0  *****************************  *.\*^\W* 

99999999999999999000900900*09****************** 

999999999990X99990999*999********************** 
9  9  9  9  t  9  0  9  •  9  ******************************  *,\W*^\* 
m9999999999t000t9*99*************************** 
999*9999999999***0***************************** 

999999999909999099***************************** 
9  9  9  9  9  9  9  9  9  0  9  . . . . 

11  19  1*  t»  1*  *1  **  /♦»#»•  1-  U  t*  *1  •»  ••  *•  ••  ^ 


FIG.  5  CELL  FLAGGING  FOR  PROBLEMS  3  AND  4 
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nant  72-ui 

m  am  mar  aatlaata.  lha  rl^  haai  sUa  af  (IS)  la  a»t  am  crxar  Saaai  siaca 
far  aaf  aatxlx  aam  ahaa  1 1  1 1  <  1  aa  feaaa 


p€LJ  <  n  M 
1  -  aaj  "  1- n  I-.  fl 


far  all  yrnSlaM  aa  haac  tte  folloalai  data: 
M  *  zadlal  dSmaamiam  of  maaM  •  56 
■  *  aertlcal  dlaaaiioa  of  andi  «  122 


•  lades  of  hoxixoacal  Idqrar  of  surface  cells 


0 

0 


•  0 


•  kooan  true  solntloa  »  -sCa,  ' s.) 

•  Initial  guess  »  0 


120 


Table  2  gives  the  associated  paraaeters  needed  tr  calculate  the  error  estl- 
aates  for  Prohleas  1*2,3  and  4.  T^les  3  and  4  give  Albrecht's  two  error  estl- 
aates  froa  equations  (10)  or  (14),  and  (11)  or  (15),  Tables  3  and  4  also  give 
the  spectral  error  estlnate  froa  equation  (16)  as  well  as  the  true  error 
1 1 4^  -  dl  1 2«  For  conparlson  all  paraaeters  are  noraallzed  by  I  1 2* 

CONCLUSION 


From  Tables  3  and  4  we  see  that  :  or  the  unevenly  soaced  nroblems,  Albrecht's 
bounds  are  over-estimates  by  more  thai  five  significant  digits,  and  that  for  the 
evenly  soaced  nroblems  they  are  over-estimates  by  about  three  significant  digits. 
This  is  due  to  the  fact  that,  as  shown  in  Table  2,  l|n  i I  2  is  large  for  the 
unevenly  spaced  ftroblems.  However,  ever  for  the  evenlv  soaced  problems,  Albrecht's 
error  bounds  are  overly  oesslmfstic  and  thus  cofioutat tonally  unusable. 
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1 

2 

3 

4 

.999949 

OMMA 

•3333WM 

.999949 

.999864 

1.980 

1.987 

1.980 

1.968 

»a,) 

.980 

.987 

.980 

.968 

9.80  Et^3 

2.53  EM 

9.80  El-3 

3.66  EM 

184. 

.354 

184. 

.354 

4.00 

4.00 

3.16 

4.00 

1.80  E46 

8.96  Bf3 

1.80  E46 

5.20  EM 

7.21  M 

3.58  EM 

5.70  EM 

5.20  EM 

6  <« 

49.0 

75.9 

49.0 

30.3 

(1)  Spectral  radlua  of  associated  Jacobi  Matrix 

(2)  Optlaue  over^'relaxatloii  factor 

(3)  Spectral  radius  of  over-relaxation  aatrlx 

(4)  Constant  troe  (10)  and  (11)  for  Albrecht's  error  estimate 

(5)  D  Is  the  diagonal  of 

(6)  Constant  from  (10)  and  (11)  where  a  ■  i)|d 

(7)  Constant  from  (10)  and  (11)  where  6  ■  oI1d^^^)!2 

(8j  5  >  p(L  )/(l-p(L  ))  for  spectral  error  estimate  (16) 

£0  (t) 

Table  2  Parameters  Needed  to  Calculate  Error  ’  stimates 
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n 

11% 

II" 

ll♦.-♦ll 

j/IUJIj 

•l!% 

(2> 

i\\2 

•IlD^^llj" 

‘lu. 

II 

»''2 

50 

4.70 

8-2 

1.38 

3.39 

E|.5 

2.01 

Ef5 

2.30 

100 

6.29 

B-3 

1.5i7 

E-1 

4.54 

E*-4 

2.92 

M 

3.08 

E-1 

150 

2.44 

B-3 

2.98 

B-2 

1.76 

M 

1.04 

Ef-4 

1.20 

E-1 

200 

8.20 

E-4 

1.92 

E-2 

5.92 

£♦■3 

2.93 

Ef3 

4.02 

E-2 

250 

2.20 

E-4 

6.80 

E-3 

1.58 

B^3 

9.68 

£*•2 

1.08 

E-2 

300 

1.38 

E-4 

4.35 

E-3 

9.99 

E^2 

5.47 

E*-2 

6.78 

E-3 

350 

5.17 

B-5 

9.88 

E>4 

3.73 

E*-2 

2.53 

B*-2 

2.53 

E-3 

400 

2.03 

E-5 

1.66 

E-4 

1.46 

B»-2 

9.24 

B^l 

9.94 

E-4 

Problea  1  —  Ihievenly  Spaced  Rectangular  Mesh 


n 

ll.JI  " 

ll♦.-♦ll 

a/ll*„llj 

»ll% 

II 

i‘'2 

-||d^ 

'‘e  1 1 

V'2 

*11% 

,112“' 

50 

2.63 

2-2 

1.78 

9.62 

Bf2 

9.59 

E+2 

2.12 

100 

7.87 

E-3 

3.58 

E-1 

2.82 

E+2 

2.81 

E+2 

6.22 

E-1 

150 

3.27 

E-3 

4.69 

E-2 

1.17 

E+2 

1.16 

E+2 

2.53 

E-1 

200 

5.93 

E-4 

5.45 

E-2 

2.12 

E+1 

2.11 

E+1 

4.68 

E-2 

250 

4.27 

E-4 

4.66 

E-2 

1.53 

E+1 

1.52 

E+1 

3.37 

E-2 

300 

4.77 

E-4 

2.63 

E-2 

1.71 

E+1 

1.70 

E+1 

3.76 

E-2 

350 

3.36 

E-4 

6.63 

E-3 

1.21 

E'fl 

1.20 

E+1 

2.66 

E-2 

400 

8.75 

E-5 

2.22 

E-3 

3.14 

3.12 

6.91 

E-3 

Problem  2  —  Evenly  Spaced  Rectangular  Mesh 


(1)  To  normalize  results  let  e  *  i)/lU  llo 

n  n  n— X  n  '  z 

(2)  Albrecht's  error  estimates  from  equations  (11)  and  (13) 

(3)  Albrecht's  error  ertimates  from  equations  (10)  and  (14) 

(4)  Spectral  error  from  equation  (16)  letting  6  =  p(I^^)/(l-plL^)) 


Table  3  Errors  for  “roblems  1  and  2 
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n 

11*. 

II  " 

1  I2 

/IhJIj 

»iis 

Trrsr 

iiT' 

50 

1.58 

E-2 

3.80 

E-1 

8.20 

E+1 

8.14 

E+1 

4.70 

E-1 

100 

3.83 

E-3 

2.42 

E-2 

1.99 

E+1 

1.98 

E+1 

1.14 

E-1 

150 

5.94 

E-4 

1.40 

E-2 

3.09 

3.06 

1.77 

E-2 

200 

1.22 

E-4 

1.24 

E-3 

6.33 

E-1 

6.29 

E-1 

3.63 

E-3 

250 

4.38 

E-5 

4.40 

E-4 

2.28 

E-1 

2.26 

d-1 

1.30 

E-3 

300 

5.56 

E-6 

3.57 

E-5 

2.89 

E-2 

2.86 

E-2 

1.65 

E-4 

350 

1.15 

E-6 

1.90 

E-5 

5.97 

E-3 

5.85 

E-3 

4.42 

E-5 

400 

2.11 

E-7 

1.76 

E-6 

1.09 

E-3 

1.09 

E-3 

6.28 

E-6 

Problem  4  —  Evenly  Spaced  Non-rectangular  Mesh 


n 

lie 

IU<« 

IIV*ll2 

/lUJU 

Bile. 

II, 

«IId^ 

/2e  1 1 

V  '2 

iiT^ 

50 

4.36 

E-2 

1.04 

2.49 

E+5 

1.90 

E+5 

2.14 

100 

5.83 

E-3 

9.40 

E-2 

3.32 

E+4 

2.68 

E+4 

2.85 

E-1 

150 

2.41 

E-3 

1.53 

E-2 

1.37 

E+4 

1.02 

E+4 

1.18 

E-1 

200 

9.03 

E-4 

1.78 

E-2 

5.16 

E+3 

3.59 

E+3 

4.43 

E-2 

250 

2.79 

E-4 

5.55 

E-3 

1.59 

E+3 

1.32 

E+3 

1.37 

E-2 

300 

1.51 

E-4 

3.53 

E-3 

8.65 

E+2 

6.38 

E+2 

7.43 

E-3 

350 

4.99 

E-5 

5.20 

E-4 

2.85 

E+2 

2.42 

e+2 

2.45 

E-3 

400 

2.15 

E-5 

1.26 

E-4 

1.23 

E+2 

9.89 

E+1 

1.05 

E-3 

Problem  3  —  Unevenly  Spaced  Non-rectangular  Mesh 


(1)  To  normalize  results  let  e^  ■  I  1  2 

(2)  Albrecht's  error  estimates  from  equations  (11)  and  (15) 

(3)  Albrecht’s  error  estimates  from  equations  (10)  and  (14) 

(4)  Spectral  error  from  equation  (16)  letting  6  ■  p(Lj^)/(l-p(Ly)) 

Table  4  Errors  for  Problems  3  and  4 
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A  surprising  computational  result  is  that  the  spectral  estimate  is  very  good 
for  all  of  the  sample  problems .  We  have  hence  decided  to  Implement  this  much 
more  practical  but  less  desirable  spectral  error  estimate  (16)  into  the  MACNOL 
computer  program. 
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